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ABSTRACT 

The  development  of  photonic  crystals  (PhCs)  has  enabled  groundbreaking  approaches  to 
mold  the  flow  of  electromagnetic  waves  with  frequency  spanning  optical  regime  to  RF  regime. 
Periodic  electromagnetic  materials  (RF  PhCs)  are  presently  one  of  the  most  rapidly  advancing 
sectors  in  the  electromagnetic  arena.  Herein,  we  demonstrated  the  applications  of  RF  PhCs  in 
improving  the  performance  of  naval  communication  system  and  decreasing  scattering  cross- 
section  of  naval  vessels.  To  this  end,  we  designed,  fabricated,  and  characterized  dielectric  devices 
based  on  the  unique  dispersion  properties  of  PhCs  to  achieve  confinement,  waveguiding, 
radiation,  and  filtering  in  the  microwave  to  millimeter-wave  portion  of  the  electromagnetic 
spectrum.  Ln  particular,  we  experimentally  demonstrated  dispersion-based  guiding  in  2D  and  3D 
PhCs.  Dispersion-based  guiding,  namely  self-collimation,  allows  for  creating  “non-diffractive 
materials”,  in  which  electromagnetic  waves  can  propagate  along  subwavelength-sized  channels 
without  divergence  while  no  specific  route  is  introduced.  In  addition,  by  combining  the  self- 
collimation  PhCs,  we  demonstrated  dielectric  grating  antennas  which  provide  stable  performance, 
compact  dimensions,  high  radiation  efficiency,  and  high  coupling  efficiency  with  other  dielectric 
waveguiding  devices.  Furthermore,  we  designed,  fabricated,  and  demonstrated  RF  channelizers 
and  filters  based  on  the  hybrid  of  properties  and  lattices  of  PhCs.  Our  work  provides  a  framework 
for  establishing  dielectric-based  naval  communication  systems,  which  are  more  secure  and  more 
efficient  compared  to  current  metallic  counterparts.  In  this  framework,  all  communication 
functions,  i.e.  generating,  sending,  transmitting,  receiving,  and  processing  RF  signals,  are  realized 
using  dielectric  devices,  which  can  be  embedded  within  structural  composite  materials  thereby 
afford  low  cross-section,  RF  functionality,  and  structural  support. 

1.0  INTRODUCTION 

Photonic  crystals  (PhCs)  were  proposed  as  ID,  2D,  or  3D  periodic,  often  dielectric  and 
artificial  structures,  which  are  designed  to  be  analogous  to  semiconductor  crystals  naturally 
formed  by  periodic  arrangements  of  atoms  or  molecules1,2.  The  propagation  of  electromagnetic 
waves  in  a  PhC  is  governed  by  Maxwell’s  equations  with  the  periodic  boundary  condition  defined 
by  the  periodic  PhC  structure,  which  are  analogous  to  Schrodinger’s  equations  for  semiconductor 
crystals.  The  solution  to  the  Maxwell’s  equations  results  in  photonic  band-structurelike  dispersion 
relationship  between  frequency,  and  wavevector,  k.  Photonic  dispersion  relationship  governs 
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the  propagation  of  photons  in  the  photonic  crystal,  mimicking  electronic  dispersion  relationship 
in  semiconductors.  Photonic  crystals  were  experimentally  demonstrated  for  effectively 
controlling  photons  in  the  same  way  as  that  semiconductor  crystals  control  electrons,  and  PhCs 
are  therefore  envisioned  as  the  semiconductors  for  photons.  The  development  of  PhCs  has 
enabled  groundbreaking  approaches  to  mold  the  flow  of  electromagnetic  waves  with  frequency 
spanning  optical  regime  to  RJF  regime. 

First  of  all,  photonic  bandgaps  (PBGs)  can  be  obtained  between  neighboring  bands  if  a 
suitable  structure  is  chosen.  For  frequencies  in  a  photonic  bandgap,  the  propagation  of  light  is 
completely  prohibited  regardless  of  any  propagation  direction 3  4  5 .  In  this  sense,  PhCs  with 
photonic  bandgaps,  namely  PBG  materials,  are  optical  insulators.  However,  if  a  point  defect  is 
introduced,  a  high-Q  cavity  for  efficient  light  confinement  and  emission  can  be  formed6,7’8. 
Design  of  cavities  responding  to  different  frequencies  can  give  rise  to  high-performance  channel- 
drop  filters9,  10 .  Similarly,  if  a  line  defect11  ’ 12  or  a  senes  of  coupled-cavities13, 14  15  16  are 
incorporated  into  such  a  PhC,  light  propagation  can  be  completely  confined  along  a  specific 
path17,18,19.  In  addition  to  nanocavities  and  waveguides,  PBG  materials  have  enabled  the  creation 
of  compact  optical  switches*"  *  and  micrometer-sized  lasers  .  Photonic  bandgap  materials  also 
found  applications  in  enhancing  efficiency  of  light  emission"13  and  microwave  antenna24. 
“Bandgap  plus  defects”  is  a  good  formula  for  solving  the  problem  of  photon  localization  and  light 
beam  routing. 

Furthermore,  PhCs  as  a  new  category  of  materials  provide  novel  optical  properties  that 
cannot  be  achieved  in  regular  materials.  In  particular,  closer  study  of  full-wavevector  dispersion 
diagrams  for  different  PhCs  has  revealed  unique  dispersion  properties,  by  which  some  first-ever 
effects  have  been  realized.  Unique  dispersion-related  phenomena  include  negative  refraction25,26, 
superprisirT  '  ,  self-collimation"  ’  ,  and  slowing  light  .  The  applications  of  these  dispersion 
properties  in  optical  regime  have  enabled  new  approaches  for  developing  high  performance 
photonic  integrated  devices.  On  the  other  hand,  by  dimensional  scaling,  the  applications  in  optical 
regime  can  find  their  corresponding  applications  in  RF  regime,  which  may  give  rise  to 
revolutionary  improvement  in  the  performance  of  RF  devices.  Photonic  crystals  allow  us  to 
engineer  control  over  the  propagation  of  electromagnetic  waves  to  an  extent  that  was  previously 
not  possible.  As  a  result,  periodic  electromagnetic  materials  (RF  PhCs)  are  presently  one  of  the 
most  rapidly  advancing  sectors  in  the  electromagnetic  arena.  Owing  to  the  tremendous  potential 
of  RF  PhC  structures,  there  is  a  plethora  of  applications  in  which  these  can  be  used. 

We  proposed  to  apply  RF  PhCs  in  improving  the  performance  of  naval  communication 
system  and  decreasing  scattering  cross-section  of  naval  vessels.  A  dominant  aspect  in  the  design 
of  most  naval  platforms  is  the  incorporation  of  functional  RF  devices  and  systems,  including 
waveguides,  routing  networks,  and  antennas.  Unfortunately,  in  many  cases  this  comes  at  the 
expense  of  reducing  structural  support  and  increasing  scattering  cross-sections.  One  means  to 
address  this  is  to  design  ship  superstructures  using  advanced  composite  materials  that  consist  of 
mainly  dielectric  materials.  In  this  case,  their  RF  scattering  properties  are  far  less  than  their  all 
metallic  counterparts.  However,  once  the  superstructure  and  masts  are  in  place,  albeit  with  low 
RF  cross-sections,  they  are  typically  outfitted  with  communication  and  radar  antennas  that  serve 
to  compromise  their  structural  integrity.  In  addition,  there  is  often  a  significant  increase  in 
scattering  cross-section  due  to  the  introduction  of  metallic  waveguides  that  serve  to  feed  and 
receive  antennas  that  are  also  mounted  on  the  superstructure.  To  overcome  this  limitation,  we 
need  to  develop  novel  routing  and  radiation  structures  that  are  largely  electromagnetically 
transparent  yet  based  on  their  unique  dispersion  properties  still  guide  and  radiate  electromagnetic 
waves. 
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In  addition  to  reducing  the  scattering  properties  of  waveguides  and  antennas  in  ship 
superstructure,  there  is  also  a  rising  problem  with  the  limitation  by  space  and  location 
considerations  due  to  the  numerous  radiating  apertures  that  are  needed.  Moreover,  antennas  are 
typically  placed  on  shipboard  locations  that  profoundly  affect  their  radiation  patterns  in 
comparison  to  their  initial  designs.  As  a  result,  antenna  performance  can  be  significantly  altered 
as  can  the  RJF  signature  of  the  ship,  due  to  the  significant  scattering  effects  of  having  a  large 
metallic  object  located  high  above  deck.  To  overcome  these  limitations,  we  propose  the  design 
and  fabrication  of  dispersion-based  antennas  that  can  be  directly  molded  into  the  dielectric 
composite  materials,  which  comprise  the  ship’s  superstructure. 

In  this  effort,  we  designed,  fabricated,  and  demonstrated  routing  and  radiating  structures  in 
the  microwave  to  millimeter-wave  portion  of  the  electromagnetic  spectrum.  These  devices  are 
based  on  the  unique  dispersion  properties  of  photonic  crystals  to  achieve  confinement, 
waveguiding,  radiation,  and  filtering.  The  main  advantage  of  this  approach  over  conventional 
techniques  is  that  dispersion  elements  are  to  a  large  extent  electromagnetically  transparent  and 
therefore  have  significantly  reduced  scattering  cross-sections.  In  addition,  such  structures  can  be 
embedded  directly  into  composite  material  structures  to  thereby  maintain  structural  integrity. 

2.0  TECHNICAL  ACCOMPLISHMENTS 


To  this  end,  we  focused  our  work  in  experimental  demonstration  of  dispersion-based  guiding 
in  2D  and  3D  photonic  crystals  and  dielectric  grating  antenna  systems.  Dispersion-based  guiding, 
namely  self-collimation,  allows  for  creating  v4non-diffractive  materials”,  in  which  electromagnetic 
waves  can  propagate  along  subwavelength-sized  channels  without  divergence  while  no  specific 
route  is  introduced.  The  applications  of  this  novel  effect  include  high-density  waveguiding  with 
low  channel  crosstalk  and  high  coupling  efficiency  with  other  dielectric  devices.  Dielectric 
grating  antennas  provide  stable  performance,  compact  dimensions,  high  radiation  efficiency,  and 
high  coupling  efficiency  with  other  dielectric  waveguiding  devices.  In  addition,  all  materials  used 
in  our  devices  are  dielectric,  which  can  be  embedded  within  structural  composite  materials 
thereby  afford  low  cross-section,  RF  functionality,  and  structural  support. 

The  research  was  carried  out  in  five  steps.  The  first  step  consisted  of  the  development  of  a 
millimeter-wave  imaging  system  for  mapping  the  field  distribution  (amplitude  and  phase)  in  a  3D 
volume.  Although  the  fabrication  of  millimeter-wave  PhCs  may  be  less  challenging  than  that  of 
near  infrared  (NCR)  or  visible  light,  the  corresponding  detection  system  for  millimeter-wave  is 
much  more  complex  than  an  NIR  or  visible  camera.  In  millimeter-wave  regime,  the  amplitude 
and  phase  of  a  coherent  wave  at  a  given  point  can  be  measured  using  a  monopole  connected  to  a 
vector  network  analyzer.  As  a  result,  the  detection  system  for  mapping  the  field  distribution  in  a 
3D  volume  was  established  by  3 -dimensionally  scanning  the  monopole.  A  particularly  attractive 
feature  of  this  system  is  its  ability  to  measure  results  in  both  amplitude  and  phase. 

The  second  step  was  to  design,  fabricate,  and  demonstrate  self-collimation  guiding  in  low 
index  contrast  photonic  crystals.  The  design  of  the  PhCs  is  based  on  their  equi-frequency 
contours  and  the  2D  finite-difference  time-domain  simulation  results.  In  the  experiments,  the 
millimeter-wave  photonic  crystals  are  fabricated  in  Rexolite  slabs  by  a  computer  numerically 
controlled  (CNC)  micro-milling  system.  Using  the  millimeter-wave  imaging  system,  we  have 
observed  and  characterized  the  self-collimation  (dispersion  guiding)  effect  for  both  the  amplitude 
and  phase  of  the  propagating  electromagnetic  wave  in  low  index  contrast  photonic  crystals. 

As  the  third  step,  we  designed,  fabricated,  and  charactenzed  a  mi  Dime  ter- wave  grating 
antenna  with  an  integrated  self-collimation  photonic  crystal.  The  photonic  crystal  provides  a 
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collimated  beam  feed  to  the  grating  antenna  by  reshaping  the  wave  front  using  its  unique 
dispersion  properties.  The  parameters  of  the  photonic  crystal  and  that  of  the  grating  are  chosen  so 
that  self-collimation  for  the  photonic  crystal  and  normal  scattering  for  the  grating  occur 
simultaneously  at  40GHz.  This  integration  constitutes  a  potential  improvement  on  the  overall 
performance  of  the  antenna.  Experimental  results  showed  good  directivity  and  high  radiation 
efficiency. 

As  the  fourth  step,  we  investigated  and  experimentally  demonstrated  self-collimation  guiding 
inside  a  3D  simple  cubic  PhC  at  microwave  frequencies.  If  the  significance  of  demonstration  of 
self-collimation  in  2D  planar  structures  lies  in  providing  a  novel  guiding  approach,  then  the 
significance  of  demonstration  of  3D  self-collimation  lies  in  opening  the  door  to  a  completely  new 
category  of  materials,  namely  “non-diffractive  materials”.  The  self-collimation  modes  were 
excited  by  a  grounded  waveguide  feeding  and  detected  by  a  scanning  monopole.  Self-collimation 
of  electromagnetic  waves  in  the  3D  photonic  crystal  was  demonstrated  by  measuring  the  3D  field 
distribution,  which  was  shown  as  a  narrow  collimated  beam  inside  the  3D  photonic  crystal,  but  a 
diverged  beam  in  the  absence  of  the  photonic  crystal. 

Furthermore,  if  sensors  and  communication  devices,  such  as  channel izers  and  correlators,  are 
incorporated  into  these  superstructures,  it  will  provide  capabilities  of  direct  capturing, 
digitalizing,  and  analyzing  electromagnetic  signals  at  the  receiving  end.  Meanwhile,  it  can  greatly 
increase  the  structural  support  and  compactness.  To  this  end,  we  designed,  fabricated,  and 
demonstrated  RE  channelizers  and  filters  based  on  the  hybnd  of  properties  and  lattices  of 
photonic  crystals.  Most  part  of  research  in  this  aspect  will  be  carried  out  in  PHASE  II. 

The  work  accomplished  in  PHASE  I  combining  that  will  be  finished  in  PHASE  II  will 
provide  a  totally  new  outlook  for  dielectric-based  RE  communication  systems  (e.g.,  naval 
platforms),  where  all  communication  functions,  i.e.  generating,  sending,  transmitting,  receiving, 
and  processing  RF  signals,  are  realized  using  dielectric  devices. 
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2.1  SELF-COLLIMATION  GUIDING  IN  A  2D  PHOTONIC  CRYSTAL 

Electromagnetic  beams  diverge  when  propagating  in  dielectric  materials  due  to  diffraction.  It 
seems  to  be  a  universal  rule  that  no  material  can  sustain  and  preserve  the  size  of  a  narrow  beam 
unless  specific  structural  routes  (e.g.  waveguides,  index  gradients32 ,  or  line-defects13)  are 
incorporated.  Recently,  development  of  PhCs1,2  has  enabled  groundbreaking  approaches  to 
control  the  flow  of  electromagnetic  waves.  Many  applications  of  PhCs  stem  from  the  existence 
of  a  photonic  bandgap  in  a  certain  frequency  range,  where  no  electromagnetic  wave  is  permitted 
to  propagate.  If  defects  are  introduced  in  PBG  materials,  localized  modes  can  be  excited.  The 
localized  modes  can  be  either  cavity  modes  for  point  defects33,34  or  guiding  modes  for  linear 
defects35.  In  this  sense,  PBG  provides  a  powerful  guiding  approach  since  the  guiding  route,  which 
depends  on  the  defects  incorporated  in  the  PhC,  can  be  extremely  narrow  and  sharply  bent33. 
However,  this  guiding  mechanism,  namely  “bandgap  plus  defects/’  also  has  drawbacks.  For 
instance,  photonic  bandgap  requires  high  index  contrast,  which  limits  the  material  selection  and 
often  results  in  large  index  mismatch  at  interfacial  boundaries.  On  the  other  hand,  recently 
substantial  interest  has  arisen  in  the  unique  dispersion  properties  provided  by  PhCs36,37'38.  A 
particular  phenomenon  for  PhCs  is  self-collimation29,30.  When  a  PhC  works  in  self-collimation 
mode,  an  electromagnetic  beam  can  propagate  within  it  without  divergence.  In  other  words, 
electromagnetic  wave  propagates  as  if  there  were  no  diffraction  in  the  material. 

The  concept  of  self-collimation  was  originally  proposed  by  Kosaka,  et  al  29,  followed  by  a 
detailed  theoretical  investigation  by  Witzens,  et  al j0  and  Chigrin,  et  al'9.  Chen,  et  al 40  and  Yu,  et 
al 41  theoretically  discussed  the  bending,  crossing,  and  splitting  properties  for  self-collimation 
PhCs.  However,  most  of  their  effort  was  focused  on  high  index  contrast  PhCs42'43.  To  decrease  the 
coupling  loss  between  common  materials  and  a  PhC,  we  focus  our  work  on  low  refractive  index 
contrast  PhCs44.  In  our  work,  we  designed,  fabricated  PhCs  and  experimentally  demonstrated,  for 
the  first  time,  self-collimation  effect  using  both  amplitude  and  phase,  for  low  index  contrast 
PhCs. 

The  dimensions  for  all  PhCs  are  scalable  with  regard  to  the  working  wavelength  because 
there  is  no  fundamental  constant  with  the  dimensions  of  length  in  Maxwell’s  equations45.  As  a 
result,  the  feature  size  for  near  infrared  (NIR)  or  visible  light  is  on  the  order  of  100’s  nanometers, 
while  that  for  millimeter-wave  (MMW)  is  on  the  order  of  100’s  microns.  Therefore,  the 
fabrication  of  MMW  PhCs  is  less  challenging  than  that  of  NIR  or  visible  light.  However,  the 
corresponding  detection  system  for  MMWs  is  much  more  complex  than  an  NIR  or  visible 
camera.  Nevertheless,  we  have  built  a  MMW  imaging  setup,  as  shown  in  Fig.  1,  based  on  an 
Agilent  85106D  network  analyzer,  which  encompasses  a  test  and  measurement  capability 
spanning  45MHz  through  1 10GHz. 

A  particularly  attractive  feature  of  this  system  is  its  ability  to  measure  results  in  both 
amplitude  and  phase.  To  this  end,  the  setup  is  composed  of  3  subsystems:  MMW  feed  system, 
motion  control  system,  and  S-parameter  measurement  system.  In  our  setup,  the  sample  is  placed 
on  the  top,  while  the  input  source,  from  port  1,  is  launched  either  through  a  monopole  antenna 
(source  1)  or  a  waveguide  (source  2),  which  works  as  a  point  source  or  windowed  beam, 
respectively.  To  avoid  the  potential  of  distributing  the  fields  during  measurement,  the  detector 
chosen  is  simply  a  1mm  monopole  antenna,  which  is  fed  back  into  port  2  of  the  network  analyzer. 
The  field  distribution  is  acquired  by  scanning  the  surface  of  the  PhCs  to  detect  the  evanescent 
waves  of  the  MMW  propagating  in  the  slab.  To  do  this,  we  assumed  the  gap  between  the  detector 
and  the  surface  is  constant,  in  which  case  the  evanescent  wave  distribution  is  proportional  to  the 
field  in  the  PhCs.  A  custom  program  was  developed  to  synchronize  the  scanning  and 
measurement.  Once  a  scanning  and  measurement  is  finished,  we  depict  the  S-parameter  value 
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with  regard  to  the  position  as  an  image.  Consequently,  each  pixel  in  the  image  corresponds  to  the 
S-parameter  (S2i  in  our  case)  value  at  a  scanning  position.  Also,  in  the  network  analyzer  the  S- 
parameters  are  given  as  complex  values,  so  we  can  obtain  both  amplitude  and  phase  distributions, 
which  is  impractical  for  NIR  or  visible  light  measurement. 


(a)  (b) 

Fig.  /.  The  MMW  imaging  setup  based  on  a  network  analyzer:  (a)  The  whole  system;  (b)  The  MMW 
feeds  and  the  detector . 

In  order  to  demonstrate  self-collimation,  a  MMW  material,  2.4mm  (or  3/32-inch)  thick 
Rexolite®  slab,  is  used  in  our  work.  Rexolite,  a  plastic-like  material  with  index  of  1.58,  has  the 
advantage  of  easy  fabrication,  low  loss,  and  low  cost.  In  particular,  it  can  be  fabricated  using  a 
CNC  micro-milling  machine.  Also,  its  low  index  provides  a  weaker  confinement  in  the  slab  for 
MMW  on  vertical  direction,  resulting  in  a  stronger  evanescent  wave  for  detection.  In  order  to 
determine  the  dielectric  constant  or  index,  we  first  launched  of  the  Rexolite  slab  a  40GHz  MMW 
and  scanned  the  Rexolite  slab  (thickness:  2.4mm)  without  any  pattern.  We  then  calculated  the 
effective  index  of  the  slab  according  to  the  guided  wavelength,  which  can  be  easily  calculated 
from  the  measured  phase  distribution.  The  effective  index,  defined  as 

nefj=ns\ndmQte,  (1) 

can  be  calculated  as 

neff=$z/$ 0  =y\,  (2) 


where  \>is  the  vacuum  wavelength  and  Xg  is  the  guiding  wavelength  along  the  slab. 

The  measured  values  are  neff TE)  =1.36,  and  n^TM)=1.21,  which  are  very  close  to  the 
calculated  values,  n^TE)  =1.36  and  ne^{TM)=1.23  based  on  the  ray-tracing  method  46.  Applying 
these  effective  indices,  we  calculated  the  EFCs  for  a  square  lattice  PhC  with  the  hole  diameter  to 
lattice  constant  ratio  2rla=0.6  using  a  plane-wave  expansion  method47.  In  Fig.  2,  the  results  of 
this  calculation  are  shown  as  the  normalized  frequency 

< j)tt=(j)a/2Trc  (3) 


and  the  normalized  wavevector 
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k„=ka/2ir. 


(4) 


where  c  is  the  speed  of  light. 


Fig .  2:  Eq id-frequency  dispersion  contours  for  2.4mm  thick  Rexolite  square  lattice  PhC(2r/a=0.6):  (a) 
for  TE  mode(ncff=  l  .36);  (b)  for  TM  mode(nefj=  1 .2 1). 

By  observing  the  EFCs,  we  anticipate  that  self-collimation  in  the  first  band  will  occur  at 
normalized  frequencies  of  to,  =0.44  for  TE,  and  a)„=0.50  for  TM  because  at  these  frequencies  the 
EFCs  are  nearly  flat  and  the  excited  mode  will  propagate  with  group  velocities  pointing  at  almost 
the  same  direction  (45°  to  the  lattice  orientation),  which  is  calculated  from48 


(5) 


In  the  2nd  band,  frequencies  shown  in  Fig.  2  cannot  excite  self-collimated  mode,  but  self- 
collimation  can  really  occur  at  higher  frequencies.  However,  parts  of  their  EFCs  are  located 
inside  of  the  so-called  “light  cone”  and  results  in  a  large  propagation  loss 49.  Light  cone  represents 
that  on  vertical  direction  of  the  slab  the  total  internal  reflection  is  just  satisfied  at  that  frequency, 
or  kn=<j)n.  Consequently,  light  cannot  be  fully  confined  at  vertical  direction  of  the  slab  due  to 
violation  of  the  total  internal  reflection  condition.  In  a  contrast,  in  our  work  we  chose  the  1st  band 
because  the  EFCs  of  the  MMW  with  these  frequencies  lie  outside  their  light  cones,  and  hence 
energy  can  propagate  with  minimum  loss. 

In  order  to  examine  the  self-collimation  effect  at  around  40.0GHz,  we  selected  the  following 
parameters:  hole  diameter  2a -2.00mm  and  lattice  constant  a=2r/0. 6=3. 33mm.  In  our  case,  two- 
dimensional  PhCs  are  fabricated  by  drilling  regularly  aligned  holes  through  the  Rexolite  slabs.  In 
the  zooming-in  inset  of  Fig.  3  (a),  the  holes  are  illustrated  in  white.  The  nearest  neighboring 
distance  between  holes  is  lattice  constant  a.  Since  the  first  band  at  the  selected  frequency  is  at  an 
angle  of  45°  to  the  lattice  orientation,  for  observational  convenience,  we  rotate  the  PhC  lattice 
orientation  at  45°  with  respect  to  the  waveguide  feed,  as  shown  in  Fig.  3  (a).  Although  the  self- 
collimation  should  be  optimal  at  ga*=0.44,  we  examine  the  frequency  o)n=0A2  to  avoid  the 
degeneracy  with  the  second  band  because,  as  shown  in  Fig.  2  (a),  eign values  only  exist  at  the  1st 
band  (solid  lines)  for  o)w=0.42,  while  eignvalues  exist  both  at  the  1st  band  and  the  2nd  band 
(dashed  lines)  for  ui„=0.44. 
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(a)  (b) 

Fig.  3:  MMW  Rexolite  PhCs  fabricated  by  the  micro-milling  method:  (a)  PhC  for  TE  self-collimation; 
(b)  PhC  for  TM  self-collimation.  Both  of  them  are  with  hole  diameter  2r=2.00mm  and  lattice  constant 
a=2r/0. 6=3. 33mm.  The  horn  has  the  dimensions  of  w6xL40xW18mm. 

According  to  the  normalized  frequency,  we  can  calculate  the  vacuum  wavelength,  Ao  —a/ 

=  7.8  mm,  and  the  corresponding  frequency  is  about  39GHz.  Based  on  the  2D  finite-difference 
time-domain  (FDTD)  method^0,  we  simulated  a  39.0  GHz  TE  MMW  from  a  monopole  source 
propagating  in  the  PhC,  as  shown  in  Fig.  4  (a).  For  TE  modes,  the  MMWs  have  three  dominant 
components,  Ex ,  Ey ,  and  Hz.  Either  the  distribution  of  magnetic  field,  Hzy  or  electrical  field, 

E  =  fiT+E;  ,  indicates  directions  of  the  power  distribution.  It  is  easier  to  measure  the  electrical 

field,  so  we  show  the  amplitude  of  the  electrical  field  in  the  simulation.  Thus,  the  simulation 
shows  the  propagation  is  primarily  in  two  orthogonal  directions,  which  indicates  the  presence  of 
self-collimation.  In  contrast,  it  is  well  known  that  the  propagation  is  isotropic  among  common 
materials. 

For  TM  modes,  the  effective  index  is  much  smaller.  As  a  result,  self-collimation  only  occurs 
over  a  very  narrow  incident  angle  range  and  the  degeneracy  with  the  second  band  is  inevitable 
because,  as  shown  in  Fig.  2  (b),  eignvalues  exist  both  at  the  first  band  (solid  lines)  and  the  2nd 
band  (dashed  lines)  for  frequencies  spanning  from  0.48  to  0.54,  whose  EFCs  are  acceptably 
straight.  We  chose  g)„=0.5,  which  corresponding  to  MMW  with  vacuum  wavelength,  Xo  =a/  u)n  = 
6.67mm,  and  frequency  45GHz  for  the  PhC  with  a=2r/0.6=3.33mm,  and  launched  a  windowed 
beam  (with  the  dimension  of  WR-22  waveguide,  w=5.6mm)  instead  of  a  monopole  as  the  input, 
as  shown  in  Fig.  3  (b).  Also,  the  input  is  expanded  through  an  integrated  horn  with  dimensions 
w6xI40x  JF18mm,  by  which  the  launched  MMW  is  expanded  to  an  18mm  windowed  beam  at  the 
distance  of  40mm.  Applying  the  2D  FDTD  method,  we  show  that,  in  Fig.  4  (b),  the  45GHz  TM 
MMW  is  self-collimated  even  for  this  PhC  with  such  a  low  index  contrast.  For  TM  modes,  the 
MMW  has  only  one  electric  component,  i.e.,  E-Ez. 

Fig.  5  shows  the  experimental  results.  From  the  images  we  can  see  that  the  Ex  distribution  is 
along  y-axis,  while  the  Ey  distribution  is  along  x-axis.  The  relative  amplitude  of  Ez  is  calculated 
from 


(6) 


E  °c 


+  E 


The  image  of  E  shows  that  the  propagation  is  primarily  along  two  directions,  namely  x-axis 
andy-axis,  and  as  a  result  the  MMW  is  self-collimated. 


x  (millimeter) 


i  jG 9 


(b) 

Fig .  4:  (a)  Simulation  of  amplitude  Hz  when  a  39GHz  TE  monopole  MMW  propagates  in  the  designed 
PhC;  (b)  Simulation  of  amplitude  E.  when  a  4 5 GHz  TM  MMW  fed  from  a  horn  propagates  in  the 
designed  PhC .  Both  PhCs  with  hole  diameter  2r-2. 00mm  and  lattice  constant  a=2r/0. 6=3. 33mm. 

To  verify  these  results  experimentally,  we  first  fabricated  (69><69+l)/2  holes  according  to  the 
simulated  structure  in  a  Rexolite  slab  and  cleaved  it  into  a  200mrnx200mm  square  with  holes 
aligned  in  the  middle  of  the  slab  as  shown  in  Fig.  3  (a).  Since  for  TE  modes  the  electric  filed 
component  is  polarized  in  the  slab,  we  bent  the  source  monopole  as  well  as  detection  monopole 
90°  to  ensure  they  were  parallel  to  the  slab  as  illustrated  in  Fig.  1(b).  We  inserted  the  source 
monopole  in  the  center  hole,  set  the  frequency  at  39GHz,  and  scanned  the  slab  with  steps 
dx-dy=  0.5mm.  To  measure  Ex ,  the  detecting  monopole  antenna  was  placed  along  x-axis,  and  Ey, 
y- axis. 

To  experimentally  validate  TM  self-collimation,  we  first  fabricated  the  PhC  structure  with 
(33*71  +  l)/2  holes  and  then  cleaved  it  into  a  slab  with  the  designed  fed-in  horn,  as  shown  in  Fig. 
3  (b).  In  our  experiment,  we  fed  this  PhC  with  a  WR-22  waveguide  with  dimensions 
5.6mmx2.8mm.  In  this  case,  the  waveguide  excites  TE0i  mode  and  electrical  field  (or  E-vector)  is 
polarized  along  the  narrower  side  of  the  waveguide.  We  rotated  the  waveguide,  making  sure  that 
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the  E-vector  was  perpendicular  to  the  slab.  The  detecting  monopole  antenna  was  straight, 
coinciding  with  the  E-vector  as  illustrated  in  Fig.  1  (b).  We  set  the  frequency  at  45GHz  and 
scanned  the  slab  with  steps  dx=dy=  0.5mm. 


(a) 


l 


(b) 


(c)  (d) 

Fig.  5:  Images  of  a  39GHz  TE  monopole  MMW  propagating  in  the  fabricated  PhC:  (a)  Image  of  Ey;  (b) 

Image  of  Ex;  (c)  Image  of  Hz  calculated  from  images  of  Ey  and  Ex  H  2  cc  ^ Je *  +  E2V  ;  (d)  Image  of  Eyfor 

39 GHz  TE  MMW  propagating  in  the  fabricated  PhC.  MMW  is  fed  from  a  waveguide  through  the 
integrated  horn . 

Fig.  6  (a,  b)  shows  the  experimental  results  for  amplitude  and  phase.  Images  of  Ez ,  both  in 
amplitude  and  phase,  show  that  the  propagation  is  along  only  one  direction  and  the  MMW  is 
indeed  self-collimated  even  for  an  effective  index  as  low  as  1.21.  Also,  due  to  the  low  contrast  the 
back-reflection  at  the  interface  to  the  PhC  is  very  small.  The  beats  shown  in  the  amplitude  might 
be  attributed  to  deviations  from  straight  lines  in  the  EFC  at  ends  and  activity  of  higher  harmonic 
inodes.  In  addition,  by  picking  up  the  narrow  central  rows,  we  can  estimate  the  propagation  loss, 
which  was  determined  to  be  O.ldB/cm. 

In  contrast,  we  fabricated  another  device  with  the  same  horn  and  dimensions,  but  without  the 
PhC  on  it.  We  repeated  the  experiment  and  obtained  the  images  for  the  amplitude  and  phase,  as 
shown  in  Fig.  6  (c,  d).  It  is  clear  that  the  propagation  in  this  device  is  not  collimated,  which  serves 
to  prove  that  the  collimation  in  our  PhC  device  is  not  originated  from  the  horn. 

Based  on  the  same  PhC  with  horn  device,  we  altered  the  source  to  39GHz  TE  MMW  and 
repeated  the  experiment.  The  result  is  shown  in  Fig.  5(d).  The  self-collimation  effect  for  TE  mode 
is  validated  again.  Using  the  same  process  as  for  TM  mode,  we  estimated  the  propagation  loss, 
which  was  determined  to  be  0.2dB/cm.  The  difference  between  the  loss  for  TE  mode  and  that  for 
TM  mode  may  be  attributed  to  the  different  feed  beam  widths:  TM  with  the  width  W0= 5.6mm  and 
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TE  with  the  width  ^=2. 8mm. 


(c)  (d) 

Fig.  6:  (a)(b)  Images  of  Ez  (amplitude  and  phase)  for  45GHz  TM  MM W  propagating  in  the  fabricated 
PhC.  MMW  is  fed  from  a  waveguide  through  the  integrated  horn /  (c)(d)  Images  of  Ez  (amplitude  and 
phase)  for  45GHz  TM  MMW  propagating  in  a  device  with  the  same  dimensions  but  without  the  PhC. 
MMW  is  fed  from  a  waveguide  through  the  horn. 

We  noticed  that  self-collimation  only  occurs  at  some  spectral  ranges.  So  we  repeated  the 
expenment  using  the  same  PhC,  but  spanned  the  frequencies  from  35GHz  to  41  GHz  for  TE,  and 
from  42GHz  to  48GHz  for  TM.  We  obtained  two  senes  of  images  showing  how  self-collimation 
begins  and  how  self-collimation  ends  with  relation  to  the  frequencies.  From  these  images,  we 
conclude  that  self-collimation  is  acceptable  in  the  ranges  of  (38.6GHz,  41.0GHz)  for  TE  mode, 
and  (43.8GHz,  46.2GHz)  for  TM  mode. 
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2.2  MILLIMETER-WAVE  DIELECTRIC  GRATING  ANTENNA  WITH 
INTEGRATION  OF  A  SELF-COLLIMATION  PHOTONIC  CRYSTAL 


Integrated  planar  gratings  have  been  widely  used  in  integrated  optics  as  distributed-feedback 
lasers,  mode  size  converting  couplers,  and  add-drop  multiplexers51,52^3,54.  Recently,  considerable 
interest  has  risen  in  the  use  of  all-dielectric  grating  antennas  due  to  their  low  loss,  low  cost,  and 
easy  integration  with  other  planar  devices55,56,57.  In  particular,  grating  antennas  for  millimeter- 
wave  applications  can  be  comparatively  compact  and  provide  acceptable  directivity.  In  addition, 
they  offer  the  capability  of  beam  scanning  and  steenng  by  varying  the  grating  parameters,  or  the 
frequency58'59.  However,  planewave-like  feedings  are  often  required  when  grating  antennas  are 
built  in;  otherwise,  their  performance  will  degrade.  Recently,  we  proposed  and  demonstrated  self- 
collimation  photonic  crystals  (PhCs),  which  can  reshape  the  wave  front  into  plane  waves60,61,62,63. 
Herein,  we  apply  this  new  technique  and  demonstrate  a  millimeter-wave  grating  antenna  with  the 
integration  of  a  self-col limation  photonic  crystal. 


Fig.  1:  The  schematic  of  the  antenna.  The  device  is  composed  of  a  square  lattice  PhC  and  a  grating 
structure.  The  parameters  of  the  PhC  (2r~2.0mm,  2r/a=0.6)  and  the  grating  (A=5.8mm,  4=2. 9mm, 
Depth=0.6mm)  are  chosen  so  that  self-collimation  for  the  PhC  and  normal  scattering  for  the  grating  occur 
simultaneously  at  40GHz.  The  thickness  of  the  Rexolite  slab  is  2.4  mm. 

Self-collimation  occurs  when  the  equi-frequency  contours  (EFCs)  are  square-like  in  shape. 
When  a  PhC  works  in  self-collimation  modes,  an  electromagnetic  beam  propagating  within  it  will 
do  so  without  divergence.  In  other  words,  an  electromagnetic  wave,  or  light,  is  channelized  as  if  it 
were  contained  within  a  linear  defect,  but  no  defect  is  present.  Furthermore,  the  size  of  the  beam 
depends  on  that  of  the  feed  size  and  can  be  large,  but  the  wave  front  of  a  self-collimation  beam 
maintains  a  flat  profile  regardless  of  the  beam  width.  This  property  provides  a  very  good 
approximation  of  a  plane  wave  and,  hence,  a  good  feed  for  a  grating  antenna.  In  contrast,  a 
standard  dielectric  waveguide  cannot  provide  a  wide  and  a  flat-phase-profile  feed  due  to  multiple 
modes. 

Fig.  7  illustrates  our  proposed  device,  which  is  composed  of  a  square  lattice  PhC  and  a 
grating  structure.  The  parameters  of  the  PhC  and  the  grating  are  chosen  so  that  self-collimation 
for  the  PhC  and  normal  scattering  for  the  grating  occur  simultaneously  for  a  given  working 
frequency.  In  the  description,  we  use  a  Cartesian  coordination  system:  the  x-axis  is  along  the 
guiding  direction,  ^y-axis  is  in  the  device  surface,  and  z-axis  is  perpendicular  to  the  device  surface, 
see  the  inset  of  Fig.  7. 
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Fig.  8:  The  equi-frequency  contour  for  the  photonic  crystal. 

The  photonic  crystal  used  in  this  work  is  similar  to  that  reported  in  last  section.  In  this  work, 
we  used  a  Rexolite  slab  with  a  thickness  of  2.4mm  and  a  bulk  refractive  index  of  1.58.  To  study 
this  planar  structure,  we  employed  the  effective  index  method,  which  simplifies  an  otherwise  3D 
problem  to  a  2D  problem  with  an  adjusted  index,  namely  the  effective  index  of  the  slab.  A  good 
approximation  for  the  effective  index  is  neff^/k0,  where  (3  and  ko  (=co/c=27rA)  are  the  propagation 
constants  in  the  slab  and  in  vacuum,  respectively.  Since  the  propagation  constants  for  TE  and  TM 
modes  are  different,  the  effective  indices  for  TE  and  TM  modes  are  also  different,  even  with 
electromagnetic  waves  of  the  same  frequency  propagating  in  the  same  slab.  For  instance,  the 
effective  indices  for  the  Rexolite  slab  are  calculated  to  be  nejf  TE)  =1.36,  and  nefj{ TM)=1.21  for  a 
propagating  wave  at  40GHz.  As  a  result,  the  dispersion  properties  of  a  photonic  crystal  for  TE 
and  TM  modes  will  also  be  different.  In  this  work,  we  only  demonstrate  the  TE  case,  where  the 
magnetic  field  is  along  z-axis.  Applying  the  effective  index  method,  we  calculated  the  equi- 
frequency  contours  of  the  dispersion  surfaces  for  2r/#=0.6  square  lattice  PhC  using  the  plane- 
wave  expansion  method,  where  a  is  the  lattice  constant  of  the  crystalline  structure,  and  2 r  is  the 
diameter  of  the  circular  holes. 

In  Fig.  8,  the  results  of  the  calculation  are  shown  for  the  normalized  frequency  tdn=fa/c  and 
the  normalized  wavevector  k,=ka/2'K.  In  the  equi-frequency  dispersion  contours,  we  anticipate 
that  self-collimation  in  the  first  band  will  occur  at  a  normalized  frequency  of  g)„=0.44  for  the  TE 
mode  because  at  this  frequency  the  equi-frequency  contour  is  nearly  straight.  In  order  to  design 
the  device  to  work  at  40  GHz,  we  choose  the  lattice  constant,  a=0A4c/f=3.3  mm  and  the  diameter 
of  the  hole,  2r=0.6a=2.0  mm. 

The  design  of  the  grating  structure  is  based  on  the  phase  matching  between  the  scattering 
wave  and  guiding  wave,  which  is  given  from  the  grating  equation64: 

„  2/r 

yx-P  =  —  m,  (7) 

A 
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where  yx  is  the  x-component  of  the  scattering  wave  vector,  (3  (=  neJy — )  is  the  propagation 

constant,  and  m  is  an  integer.  See  Fig.  9(a).  To  achieve  normal  scattering,  7=0.  Thus,  one 
straightforward  choice  is  that  m=- 1,  resulting  A  =  }Jnejj=  5.6mm  for  a  TE  mode  in  our  slab 
operating  at /=40GHz.  Appling  this  grating  period  in  a  2D  finite-difference  time-domain  (FDTD) 
simulation  and  carefully  tuning  the  parameters,  we  verified  and  optimized  the  design  by 
evaluating  the  effect  of  the  grating  on  nejy.  In  particular,  the  grating  period  should  be  tuned  up 
because  the  effective  thickness  of  the  slab  decreases  and,  hence,  the  effective  index  decreases 
when  the  grating  pattern  is  imposed  on  a  blank  slab.  The  simulation  results  are  shown  in  Fig. 
9(b).  Accordingly,  we  selected  another  structure  with  the  following  excitation  parameters:  slab 
thickness  T=2A  mm,  grating  period  A=5.8  mm,  grating  period  aspect  ratio  A/A=0.5,  grating 
depth  d- 0.6  mm,  source  wavelength  \=7.5  mm  (f=40  GHz),  and  a  source  polarization  of  TE 
mode. 
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Fig.  9:  (a)  The  illustration  showing  the  principle  of  the  antenna  grating .  (b)  The  phase  distribution  by  the 
two-dimensional  FDTD  simulation. 

In  our  experiment,  we  fabricated  the  device  using  a  computer  numerically  controlled  (CNC) 
router.  A  tapered  planar  structure  is  incorporated  to  efficiently  feed  a  launched  millimeter-wave 
into  the  device.  We  employed  our  millimeter-wave  imaging  system61,65,66,67  as  shown  in  Fig.  10  to 
map  the  near  field  distribution  three-dimensionally.  This  system  is  mainly  composed  of  an 
Agilent  85 106D  vector  network  analyzer  and  an  XYZ  scanner.  In  our  setup,  the  sample  is  placed 
on  the  platform,  while  the  input  source,  from  port  1,  is  launched  through  a  WR-22  waveguide.  To 
avoid  the  potential  of  distributing  the  fields  during  measurement,  the  detector  chosen  is  simply  a 
1mm  monopole  antenna,  which  is  fed  back  into  port  2  of  the  network  analyzer.  The  field 
distribution  is  acquired  by  scanning  along  the  surface  of  the  PhCs  to  detect  the  propagation  and 
evanescent  waves  of  the  MMW  propagating  in  the  slab.  A  custom  program  was  developed  to 
synchronize  the  scanning  with  field  measurement.  Once  the  scanning  and  measurement  is 
finished,  we  depict  the  S-parameter  value,  S2\  in  our  case,  with  regard  to  the  position  as  an  image. 
Consequently,  each  pixel  in  the  image  corresponds  to  the  S-parameter  value  at  each  scanning 
position. 


14 


Fig.  10:  Experimental  setup  for  charactering  the  device. 


To  visualize  the  field  better,  we  kept  the  value  of  the_y-axis  constant  and  made  a  scan  along 
the  vertical  (x,  z )  plane.  Fig.  1 1  shows  the  measured  amplitude  and  phase  distributions.  From  the 
phase  distribution,  we  can  see  that  the  scattering  waves  are  along  the  z-axis,  which  is  the  normal 
of  the  device  surface  and  in  good  agreement  with  the  simulation  results  shown  in  Fig.  9  (b). 
However,  the  amplitude  distribution  shown  is  not  yet  uniform.  This  issue  may  be  resolved  by 
carefully  chirping  the  aspect  ratio  of  the  grating.  Such  chirped  gratings  are  still  under 
development  and  will  be  reported  in  future  publications. 
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Fig.  !  1 :  Field  distribution  (phase  and  amplitude)  along  (xz)  plane. 

In  addition,  we  also  kept  value  of  the  z-axis  a  constant  (where  the  distance  to  the  surface  is 
constant)  and  made  several  scans  along  the  horizontal  (x,  y)  plane.  Fig.  12(a)  shows  the  amplitude 
and  phase  distribution  along  the  horizontal  (x,  y)  plane,  when  the  detector  is  7  mm  over  the  device 
surface.  Since  the  distance  is  quite  large  compared  with  the  wavelength,  most  of  the  waves  sensed 
by  the  detector  are  scattered  waves,  which  propagate  into  far  field.  This  emission  of  power  makes 
the  grating  an  antenna  (The  slight  phase  variation  along  x-axis  may  be  attributed  to  fabrication 
errors  in  the  grating  structure.)  On  the  other  hand,  Fig.  12(b)  shows  the  amplitude  and  phase 
distribution  along  the  horizontal  (x,  y)  plane  when  the  detector  is  2mm  over  the  device  surface. 
Since  the  distance  is  small  compared  with  the  wavelength,  most  of  the  waves  sensed  by  the 
detector  are  evanescent  waves,  which  exist  in  the  device  and  propagate  as  guided  waves.  The 
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phase  distribution  shown  in  Fig.  12(b)  is  a  series  of  parallel  lines.  This  is  attributed  to  the 
realization  of  the  self-collimation  photonic  crystal.  The  amplitude  distribution  shown  in  Fig.  12(a, 
b)  decays  because  the  wave  is  leaking/scattering  through  the  grating. 
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Fig.  12:  Field  distribution  along  the  surface  of  the  structure  (xy  plane),  (a)  The  detector  is  7mm  over  the 
surface .  (b)  The  detector  is  2mm  over  the  surface. 
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2.3  SELF-COLLIMATION  IN  A  3D  PHOTONIC  CRYSTAL 


The  demonstration  of  self-collimation  in  2D  PhCs38'29’68,  provides  a  significant  guiding 
approach,  where  nevertheless  the  vertical  confinement  of  light  is  still  achieved  by  index  guidance. 
Numerical  investigations  of  self-collimation  in  3D  PhCs  were  firstly  reported  in  Refs.  69  and  70, 
and  recently  in  Refs.  71  and  72,  offering  tantalizing  new  evidence  that  self-collimation  in  3D 
PhCs  should  be  possible,  and  opening  the  door  to  producing  a  totally  new  category  of  materials, 
namely  “non-diffractive  materials”.  However,  no  experimental  work  has  been  reported  to  date.  If 
the  significance  of  demonstration  of  self-collimation  in  2D  planar  structures  lies  in  providing  a 
novel  guiding  approach,  then  the  significance  of  demonstration  of  3D  self-collimation  lies  in 
opening  the  door  to  a  completely  new  category  of  materials,  namely  “non-diffractive  materials”. 
The  challenges  are  mostly  attributable  to  the  difficulty  of  3D  structure  fabrication,  and  3D  field 
distribution  measurement.  As  far  as  we  know,  there  has  never  been  experimental  work  reported  to 
measure  field  distribution  inside  a  3D  PhC  at  all.  In  this  work,  we  have  overcome  these 
challenges.  Herein,  we  demonstrate  a  3D  PhC  as  a  new  material,  in  which  light  can  propagate 
without  divergence  while  no  specific  route  is  introduced.  A  3D  PhC  with  desired  dispersion  was 
fabricated,  and  self-collimation  was  experimentally  venfied  by  mapping  the  field  distribution 
inside  the  volume  of  the  3D  PhC. 


Fig.  13:  (a)  The  simple  cubic  PhC  fabricated  by  a  high  precision  computer-controlled  machine.  It  has 
1 8><l  5  *36units.  (b)  The  photonic  band  diagram  for  the  PhC. 


Fig.  13(a)  shows  the  3D  PhC  we  employed  in  this  work.  It  is  a  simple  cubic  lattice  with  the 
unit  cell  shown  in  the  inset.  The  hole  diameter  is  2/ -3.2mm,  the  lattice  constant  is  a=5.0mm  and 
the  microwave  material  has  a  dielectric  constant  er=30  (HiK500F,  obtained  from  Emerson  & 
Cumming).  The  overall  dimensions  are  IT85mmx//76mmxl  170mm.  The  propagation  of 
electromagnetic  waves  in  a  PhC  is  governed  by  Maxwell’s  equations  with  the  periodic  boundary 
condition  defined  by  the  PhC  structure.  In  designing  the  structure,  we  used  a  plane-wave 
expansion  method3,73  to  solve  the  Maxwell’s  equations  and  obtain  the  eigenvalues  of  frequency / 
for  a  given  wavevector  k.  Fig.  13(b)  shows  the  /- k  diagram  around  the  perimeter  of  the 
irreducible  first  Brillouin  zone,  also  called  the  photonic  band  diagram.  According  to  the  diagram, 
a  small  bandgap  appears  between  14.4GHz  and  15.0GHz,  but  that  is  not  our  center  of  attention. 
Instead,  we  focus  on  a  frequency  region  centered  at  12.3GHz,  where  self-collimation  is  expected. 
Fig.  14(a)  shows  the  equi-frequency  surfaces  (EFSs)  30,74  at  12.3GHz,  which  connect  the 
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corresponding  wavevector  k  values  with  eigenvalue  /=12.3GHz.  Notice  that  the  EFSs  are 
depicted  in  k-space  and  two  degenerate  bands  exist.  The  third  and  fourth  bands  are  shown  in  red 
and  blue,  respectively  (The  third  band  is  also  depicted  separately  in  Fig.  14(b)).  Since  group 
velocity  (and  hence  the  direction  of  power  flow)  is  calculated  by  the  gradient  of  frequency  in  k- 
space,  i.e.,  =  2/zVA/(k),  power  propagation  direction  is  perpendicular  to  the  EFSs.  Although 

two  bands  are  degenerate,  both  EFSs  are  nearly  flat  and,  hence,  3D  self-collimation  modes  are 
supported. 


Fig.  14:  (a)  The  equi-frequency  surfaces  for  frequency  f=  1 2.3GHz.  The  third  and 
fourth  bands  are  degenerate .  (b)  The  equi-frequency  surface  of  the  fourth  band,  (c)  The 
equi-frequency  surfaces  shown  as  a  series  of  contours  for  different  L  components.  The 
black  circle  on  each  diagram  is  the  dispersion  contour  for  air . 
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A  high  precision  computer  controlled  machine  was  employed  to  fabricate  devices  for 
microwave  and  millimeter-wave  applications25,26,67,68.  This  3D  simple  cubic  PhC  was  fabricated 
by  drilling  aligned  holes  along  the  three  orthogonal  axes  in  six  85mmx85mmx25mm  pieces  of 
the  dielectric  material.  These  pieces  were  aligned  and  stacked  to  form  the  complete  PhC.  The 
fabrication  precision  was  controlled  to  ±0.1  mm,  which  is  very  small  compared  with  the  lattice 
constant  and  the  wavelength  of  electromagnetic  waves  we  will  use.  To  perform  the  experiment, 
we  established  a  3D  microwave  imaging  system67.  The  detector,  a  monopole  1.8mm  in  diameter 
with  a  1-mm  length  of  core  exposed  at  the  end,  was  mounted  on  an  XYZ  scanner.  The  fixed 
source  was  a  waveguide  with  a  flat  copper  ground  plane  at  the  end,  and  electric  component 
polarized  along  the  vertical  axis.  This  configuration  of  the  source  is  equivalent  to  a  magnetic 
dipole  source75,  which  can  excite  the  self-collimation  modes  in  the  third  band.  The  source  and  the 
detector  were  connected  to  a  vector  network  analyzer,  through  which  the  electric  field  distribution 
can  be  measured  at  any  given  point.  A  custom  program  was  developed  to  synchronize  the  motion 
of  the  detector  and  the  measurement  of  the  corresponding  field  distribution.  The  diameter  of  the 
detector  (1.8mm)  is  smaller  than  that  of  the  holes  (3.2mm).  So  we  can  insert  the  detector  into  the 
holes  to  map  the  field  distribution  inside  all  holes  of  the  PhC  following  a  point-by-point  and  hole- 
by-hole  process.  The  presence  of  the  detector  introduces  a  defect,  which  will  disturb  the  guiding 
route  defined  by  linear  boundaries,  e.g.  defect-guiding.  However,  in  this  case  the  guiding  route  is 
defined  by  the  widely  spreading  periodicity  of  the  PhC.  Its  influence  on  the  field  distribution  is 
negligible.  This  is  also  verified  by  the  experimental  results  presented  in  the  following  text. 

Ln  our  experiment,  we  varied  the  frequency  from  11.0GHz  to  17.0GHz  at  intervals  of 
0.1  GHz,  and  measured  the  electric  field  at  75  points  along  the  length  of  each  hole  separated  by 
1.0mm  steps.  As  a  result,  the  electric  field  amplitude  distribution  was  obtained  as  a  3D  matrix 
with  X75xY17xZ34  elements76  for  each  frequency.  To  estimate  the  beam  size  within  the  PhC,  we 
used  a  standard  definition  of  full  width  at  half  maximum  (FWHM),  where  we  reduced  the  effect 
of  random  variation  by  averaging  the  five  highest  values  in  each  xy  slice  as  the  reference 
maximum  value  of  that  slice.  To  show  the  change  of  the  beam  size  along  the  propagation 
direction,  we  determined  the  half-maximum  of  the  intensity  in  each  xy  slice  along  the  z-axis  and 
rendered  all  of  the  half-maximum  intensity  contours  as  a  3D  surface.  As  a  result,  the  beam  profile 
along  the  length  of  the  PhC  is  indicated  by  the  enclosed  part  of  the  3D  surface.  Fig.  15(a)  shows 
the  beam  profile  at/=12.2GHz,  which  is  a  narrow,  collimated  beam.  The  intensity  distributions  of 
slices  located  at  z=15mm  and  z=165mm  are  also  shown  in  the  insets  of  Fig.  15(a),  respectively. 

We  calculated  the  beam  diameters  on  all  slices  along  the  z-axis,  and  depicted  them  in  Fig. 
15(b).  Curve  fitting  shows  that  the  diameter  of  the  beam  increases  2mm  (from  8.5mm  to  10.5mm) 
at  the  propagation  distance  of  165mm.  Close  examination  reveals  the  increase  is  mostly  due  to  a 
jump  at  z=80mm,  where  considerable  lattice  misalignment  exists.  Optimizing  the  fabrication 
accuracy  will  improve  the  self-collimation.  In  contrast,  at  the  same  propagation  distance  a 
Gaussian  beam77  (with  the  waist  size=8.5mm)  theoretically  expands  its  waist  into  150mm  in  the 
air  and  29mm  in  a  material  with  er=30,  which  are  shown  in  Fig.  15(b)  as  the  red  and  blue  curves, 
respectively 

Detailed  analysis  demonstrated  that  this  self-collimation  phenomenon  was  observed  in  the 
frequency  range  12. lGHz~12.9GHz,  consistent  with  the  flat  EFSs  calculated  for  the  third 
dispersion  band  at  these  frequencies.  For  comparison,  we  repeated  the  same  measurement  and 
analysis  for  the  electromagnetic  wave  propagation  in  the  air  by  removing  the  PhC  while  keeping 
all  other  instruments  and  their  states.  Supplement  3  shows  a  divergent  beam  profile  in  the  air 
when  the  source  is  set  at/=12.2GHz.  The  discrepancy  of  the  profiles  from  ideal  divergent  beams 
is  due  to  the  reflection  from  environments.  No  collimated  beam  in  the  air  is  observed  at  any 
frequency  in  the  range  we  measured.  Therefore,  the  comparison  of  Fig.  15(a)  and  Supplement  3 
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further  confirms  that  it  is  the  presence  of  3D  PhC  that  results  in  the  collimation  of 
electromagnetic  beam.  In  other  words,  the  PhC  cancels  the  diffraction  effect  and  functions  as  a 
“non-diffractive  material”. 


Fig.  15:  (a)  The  FWHM  beam  profile  in  the  PhC  at  f—  1 2.2GHz.  (b)  The 
evolution  of  beam  diameter  in  the  PhC  and  common  materials. 


To  observe  how  power  is  transported  inside  the  PhC,  we  integrated  the  intensity  over  the  first 
xy  slice  and  used  it  as  input  power,  and  the  last  slice  as  output  power.  The  propagation  attenuation 
in  the  PhC  can  be  calculated  as  the  ratio  between  output  and  input.  For  different  frequencies,  we 
obtained  different  attenuations.  Fig.  16(a)  shows  the  relationship  between  the  attenuation  and 
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frequency.  From  this  diagram,  one  can  see  that  the  largest  attenuation  (79dB)  occurs  inside 
bandgap  (15.1GHz)  while  the  smallest  attenuation  (15dB)  occurs  when  the  beam  is  self- 
collimated  (12.4GHz).  Green  and  red  shaded  bars  indicate  the  minimum  and  maximum 
attenuation  regions,  respectively,  which  also  coincide  the  self-collimation  and  bandgap  ranges, 
respectively.  In  bandgap  ranges,  most  of  power  dies  out  in  the  first  several  lattices.  In  the  self- 
collimation,  the  power  (in  dB)  decays  approximately  linearly  along  the  z-axis,  as  shown  in  Fig. 
16(b),  except  for  considerable  insertion  loss  m  the  coupling  process.  Similar  issue  also  occurs  in 
2D  self-collimation  which  was  reported  in  Ref.  78.  This  may  be  solved  by  index  matching  or 
other  novel  approaches. 


Fig.  16:  (a)  The  relationship  between  the  attenuation  and  working  frequency,  (b)  The  attenuation  along  the 
2-axis  at  f=  1 2. 4 GHz. 
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2.4  HYBRID  LATTICE  PHOTONIC  CRYSTAL  CHANNELIZERS 


Among  the  properties  of  photonic  crystals,  the  unique  dispersion  and  the  bandgap  are  two 
representative  characteristics  that  allow  for  designing  devices  by  engineering  the  dispersion 
property  or  introducing  point-  and  line-defects  into  PhCs.  The  dispersion  property  offers  the 
electromagnetically  transparent  property  and  thereby  significantly  reduced  the  scattering  cross- 
sections  of  devices.  On  the  other  hand,  the  bandgap  property  provides  various  functional  devices, 
including  channelizers,  correlators,  and  filters.  Generally,  these  devices  are  based  on  a  single 
lattice  or  single  property,  such  as  line-defect  waveguides  based  on  a  triangular  lattice,  and  self- 
collimation  based  on  a  rectangular  lattice.  However,  for  the  goal  of  eventually  incorporating  of 
different  photonic  crystal  based  devices  into  subsystems  or  systems,  hybrid  of  different  lattices  is 
an  essential  task.  In  particular,  it  enables  one  to  design  more  complex  devices  by  combining 
different  properties  of  the  PhCs.  Therefore,  we  proposed  and  demonstrated  devices  based  on  the 
combination  of  these  two  approaches  that  allow  for  both  low  scattering  cross-section  and 
functional  devices. 

The  proposed  hybrid  lattice  structure  is  based  on  the  hybrid  of  the  rectangular  and  triangular 
lattices  of  air  holes  on  a  slab,  as  shown  in  Fig.  17.  The  second  pass  band  of  the  rectangular  lattice 
is  designed  to  be  in  the  bandgap  of  the  triangular  lattice,  as  shown  in  Fig.  18. 
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Fig.  17:  The  hybrid  lattice  structure  with  a  rectangular  lattice  of  air  holes  in  the  middle  surrounded  by  a 
triangular  lattice. 


In  this  case,  the  triangular  lattice  is  used  both  as  an  efficient  boundary  to  confine  light  in  the 
in-plane  direction  and  as  a  base  to  introduce  other  devices  such  as  cavities  formed  by  point 
defects.  On  the  other  hand,  the  wave  propagating  in  the  rectangular  lattice  can  be  manipulated  by 
designing  its  dispersion  properties.  Therefore,  both  the  bandgap  and  the  special  dispersion 
property  of  the  PhC  are  utilized  in  the  structure.  One  straightforward  application  of  this  hybrid 
lattice  structure  is  the  low  cross-section  waveguide. 
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Fig.  18  :  Dispersion  diagrams  of  the  triangular  and  rectangular  lattices 


To  demonstrate  the  hybrid  lattice  waveguide,  we  first  introduce  the  unique  dispersion 
property  of  the  photonic  crystal.  The  dispersion  surfaces  of  photonic  crystals  gives  one  the  ability 
to  predict  the  direction  of  light  propagating  in  a  PhC  by  studying  the  shape  of  its  equi-frequency 
surfaces  (EFCs).  Unlike  ordinary  materials  which  have  ellipsoidal-shaped  EFCs  as  shown  in 
Fig.  19(a),  PhCs  can  exhibit  a  wide  variety  of  EFC  shapes,  which  can  be  modified  by 
manipulating  the  parameters  of  the  PhC  lattice,  such  as  the  lattice  type,  pitch,  and  fill-factor.  If 
the  EFC  is  a  square-like  shape,  shown  as  the  dashed  line  in  Fig.  19(b),  the  directions  of  the  group 
velocity  are  mostly  limited  to  the  two  directions  perpendicular  to  the  square  edges.  If  light  is 
launched  toward  one  edge  of  the  EFC  in  a  wide  range  of  angles,  it  is  only  able  to  propagate  in  a 
narrow  range  of  angles  in  the  material.  In  this  situation,  light  propagation  in  the  PhC  is  self- 
collimated. 


Fig.  19:  EFCs  of  the  rectangular  lattice:  (a)  first  hand  and  (b)  second  band. 


If  EFCs  were  a  perfect  square,  light  would  propagate  in  the  photonic  crystal  without  any 
divergence.  However,  EFCs  have  curved  edges  for  real  structures.  As  a  result,  the  beam  always 
diverges  to  some  extent  as  it  propagating  in  the  lattice,  as  shown  in  Fig.  20.  To  deal  with  this 
problem,  we  created  a  hybnd  lattice  structure  with  a  triangular  lattice  as  the  boundary  to  confine 
the  light  propagating  in  the  rectangular  lattice  area.  In  this  configuration,  if  the  wavelength  of 
light  is  in  the  bandgap  of  the  tnangular  lattice,  the  divergent  light,  resulting  from  the  curved  edges 
of  the  EFCs,  is  reflected  back  to  the  waveguide.  Therefore,  the  wave  propagating  in  the 


23 


rectangular  lattice  is  highly  guided  without  divergence.  In  order  to  illustrate  the  efficiency  of  the 
hybrid  structure,  TE  polarized  waves  enter  a  rectangular  lattice  and  a  hybrid  lattice  waveguide 
with  the  same  length  simulated  using  the  2D-FDTD  method,  as  shown  in  Fig.  21 . 
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Fig.  20:  The  2D-FDTD  simulation  of  light  propagating  in  the  PhC  with  a  square-like  EFC. 


Fig.  21 :  The  2D-FDTD  simulation  of  light  propagating  in  the  PhC  with  a  square-like  EFC. 


For  the  hybnd  lattice  waveguide,  the  rectangular  lattice  with  three  periods  in  width  is  used  as 
the  waveguide.  On  each  side  of  it  is  the  triangular  lattice  serving  as  the  boundary.  In  comparison 
to  the  single  rectangular  lattice,  the  hybrid  lattice  waveguide  reveals  better  confinement  by 
increasing  the  transmission  of  20%  within  the  propagation  distance  of  600  periods.  This 
simulation  results  indicated  the  advantage  of  the  hybrid  lattice  structure  over  single  lattice 
structure  by  further  increase  the  transmission.  Moreover,  if  incorporating  point  defects  into  the 
triangular  lattice,  we  can  design  wavelength  selective  devices,  such  as  WDMs  and  channelizers. 

The  photonic  crystal  WDM  in  telecommunication  wavelength  has  been  demonstrated  to 
reduce  the  device  size  hundreds  and  thousands  of  times  more  than  currently  used  WDMs.  In 
addition  to  this  compactness,  the  motivation  for  using  this  type  of  demultiplexer  is  the  possibility 
of  integration  with  other  PC  devices  on  a  single  substrate,  which  makes  them  a  potential 
candidate  for  communication  systems,  spectroscopy,  and  sensing.  However,  the  high  loss  of  the 
line-defect  channel  waveguide  precludes  its  practical  applications.  Herein,  we  propose  an 
alternative  hybrid  lattice  PhC  WDM  device  with  a  low-loss  semi-self-collimated  waveguide  as 
the  channel  waveguide.  In  particular,  the  dropping  characteristics  are  controlled  by  manipulating 
both  the  dispersion  properties  of  the  rectangular  lattices  and  the  parameters  of  the  cavities. 
Therefore,  this  configuration  is  more  favorable  than  general  line-defect  based  WDMs,  in  which 
the  dropping  characteristics  can  only  be  controlled  by  manipulating  the  parameters  of  the 
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triangular  lattice.  To  form  the  dropping  channels,  point  defects  in  the  vicinity  of  the  channel 
waveguide  are  introduced  in  the  triangular  lattice  region.  In  this  configuration,  if  waves  propagate 
in  the  channel  with  the  frequencies  coinciding  with  the  resonance  frequencies  of  the  cavities,  it 
might  drop  to  the  cavities,  whereas  the  rest  of  the  frequencies  continue  to  propagate  in  the 
channel  waveguide. 

As  in  telecommunication  applications,  channelizers  are  important  devices  in  millimeter-wave 
communication  systems  as  well.  Therefore,  this  design  has  been  scaled  to  millimeter-wavelength 
in  an  attempt  to  design  high-compact  millimeter-wave  communication  devices.  The  hybrid  lattice 
MMW  channelizer  is  composed  of  two  rows  of  rectangular  lattice  as  the  channel  waveguide  with 
the  triangular  lattices  as  the  boundary.  Both  the  rectangular  and  the  triangular  lattice  are  designed 
with  the  same  lattice  constant:  a  =9.7  mm  and  r/a  =  0.35 ,  as  shown  in  Fig.  23.  Point  defects  are 
introduced  to  the  triangular  lattice  as  the  dropping  channels.  The  air-hole  sizes  around  the  cavities 
are  adjusted  to  make  the  desired  frequencies  drop  into  different  cavities. 


Fig.  22:  The  schematic  picture  of  the  designed  hybrid  lattice  WDM,  where  a  rectangular  lattice  serves 
as  the  channel  waveguide  and  triangular  lattices  serve  as  the  boundary.  Point  defects  are  introduced 
in  the  triangular  lattice  region  as  the  cavities. 


To  experimentally  validate  it,  the  hybrid  lattice  PhC  channelier  was  fabricated  on  a  2.4  mm 
thick  HK  12  slab  with  a  dielectric  constant  of  12.5  using  a  computer  numerically  controlled 
(CNC)  machine.  To  simplify  the  structure  for  initial  testing,  the  fabricated  device  only  consisted 
of  one  dropping  channel.  The  fabricated  channelizer  was  characterized  by  using  an  Agilent 
85106D  network  analyzer  based  system.  A  TE-polarized  MMW  with  the  wavelength  swept  from 
10  GHz  to  14  GHz  was  fed  into  the  device.  A  monopole  detector  was  used  to  map  the  surface  of 
the  device  to  detect  the  surface  scattered  field.  The  detected  signals  were  feedback  to  the  network 
analyzer  and  processed  for  the  resulting  of  the  surface  scattered  field  pattern  of  the  device.  To 
verify  the  performance  of  the  device,  the  measured  results  were  compared  with  the  2D-FDTD 
simulation  results,  which  are  consistent  with  one  another,  as  shown  in  Fig.  23.  In  both  cases,  the 
results  show  that  with  the  10.8  GHz  wave  entering,  the  wave  drops  to  the  cavity,  whereas,  with 
the  12  GHz  wave  entering,  the  wave  propagates  along  the  waveguide. 
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Fig.  23:  (a)  and  (c),  experimental  measured  surface-scattered  fields  at  12  and  10.8  GHz;  (b)  and  (d), 
2D-FDTD  simulations  of  steady-state  Hz  fields  at  12  and  10.8  GHz;  (e),  fabricated  HLPWDM ;  (f),  2D- 
FDTD  calculated  transmission  spectmm  at  the  cavity . 


To  further  study  the  mechanism  of  this  hybrid  lattice  structure,  we  calculated  the  dispersion 
diagram,  field,  and  phase  distribution  of  this  structure  using  the  2D-PWM  method,  shown  in 
Fig. 24.  Multimodes  show  up  in  the  bandgap,  where  10.8  GHz  frequency  is  belonged  to  a  flat 
band,  or  a  cavity  mode.  Meanwhile,  the  12  GHz  belongs  to  a  guiding  mode  of  the  structure.  The 
2D-PWM  simulation  results  are  consistent  with  both  the  experimental  and  2D-FDTD  simulation 
results  as  shown  in  Fig.  23. 

Encouraged  by  the  promising  experimental  results,  we  designed  a  compact  1  x  4  channelizer 
with  a  size  of  1 0A,  using  the  2D-FDTD  method.  The  structure  consists  of  a  channel  waveguide 
composed  of  two  rows  of  rectangular  lattice  and  drop  channels  formed  by  coupled  cavities  with 
different  cavity  sizes,  as  shown  in  Fig.  25(a).  The  coupled  cavity  channels  are  formed  by  the 
adjacent  point  defects  to  increase  the  dropping  efficiency  as  well  as  to  guide  the  dropped  waves. 

The  resonant  frequency  of  each  channel  is  adjusted  by  increasing  the  radius  of  the  air  holes 
around  the  cavities  as  well  as  introducing  reduced-size  air  holes  into  the  cavities.  Some  of  the 
radii  of  the  air  holes  around  the  defects  are  enlarged  to  8.1,  8.8,  9.1  and  8.2  mm  respectively  to 
ensure  the  resonant  wavelengths  for  the  respective  channels  are  as  10.8,  11.1,  1 1.4  and  12.2  GHz, 
respectively.  Fig.  25(b)  shows  the  spectra  detected  at  each  branch  channel  with  the  detectors 
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placed  in  them.  The  sharp  peaks  presented  at  the  spectra  indicate  these  wavelengths  are  dropped 
into  their  respective  branch  channels  from  the  channel  waveguide. 


Fig.  24:  2D-PWM  calculation  of  the  dispersion  diagram  of  the  hybrid  lattice,  structure  field  patterns, 
and  phase  at  the  10.8  and  the  12  GHz. 


Frequency  (GHz) 

Fig.  25:  Schematic  of  the  1x4  WDM  based  on  the  hybrid  lattice. 


To  validate  the  design,  the  structure  is  simulated  with  a  TE-polarized  wave  coupling  into  the 
rectangular  lattice  channel  waveguide.  The  steady  state  results  are  consistent  with  the  spectra. 
The  desired  wavelengths  are  dropped  to  their  respective  channels,  shown  in  Fig.  26.  The 
simulation  and  experimental  results  clearly  demonstrated  the  ability  of  the  proposed  device  to 
drop  different  frequencies  from  a  channel  waveguide  to  multiple  channels.  Based  on  the  same 
mechanism,  multichannel  demultiplexers  can  be  realized  with  careful  design  of  the  dispersion 
surface  of  the  rectangular  lattice  and  the  parameters  of  the  cavities. 
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Using  a  rectangular  lattice  as  the  channel  waveguide  for  this  hybrid  lattice  WDM  has  several 
advantages  in  addition  to  the  ultra-small  size.  First,  it  gives  more  freedom  for  designing  more 
efficient  and  complex  functional  devices  by  engineering  both  the  point  defects  and  the  dispersion 
surfaces  of  the  lattices.  In  addition,  the  similarity  between  the  hybrid  structures  and  the  CROW 
structures  proposed  by  S.  Fan’s  group  for  light-stopping  devices  may  suggest  a  new  way  to 
realize  the  phenomenon.  Interestingly,  if  the  wavelength  is  within  both  the  bandgaps  of  the 
triangular  and  rectangular  lattices,  light  can  be  trapped  in  the  cavity  resulting  in  high  Q-value  of 
the  cavity  and  a  long  lifetime  for  the  defect  modes.  Finally,  since  coupling  is  very  sensitive  to 
both  the  lattice  and  cavity  parameters,  it  is  very  easy  to  control  the  adding/dropping  of  desired 
frequencies  by  incorporating  nonlinear  materials  into  the  structure. 
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Fig.  26:  2D-FDTD  calculated  steady  state  Hz  fields  of  I  x  4  hybrid  lattice  PhC  WDM  with  the  10.8 
GHz ,  1 1.  /  GHz,  1 1 .4  GHz,  and  12.2  GHz  wavelengths  of  light  dropped  at  different  channels. 


As  demonstrated,  PhCs  provide  a  great  opportunity  to  implement  various  wavelength 
selective  devices,  including  channelizers.  These  devices  are  based  on  2D  photonic  crystals, 
wherein  the  waves  enter  the  slab  in-plane.  However,  with  the  wave  entering  the  slab  out-of-plane, 
the  transmission  spectrum  exhibits  guided  resonance  phenomenon  that  is  significantly  different 
from  the  well-studied  in-plane  incident  cases.  Herein,  we  present  a  displacement  tunable  filter 
based  on  this  guided  resonance  phenomenon. 

Most  studies  of  2D-PhCs  are  focused  on  light  that  is  confined  in  PhC  slabs  or  propagating  in¬ 
plane  of  the  slabs.  Relatively  little  research  has  been  devoted  to  studying  light  normally  incident 
to  a  PhC  slab,  in  which  the  observed  phenomena  are  significantly  different  from  the  well-studied 
in-plane/parallel  incident  cases.  A  typical  feature  of  a  case  in  which  light  is  normally  incident  is 
the  transmission  spectrum  characterized  by  the  sharp  “Fano-like”  resonance  feature,  or  the  guided 
resonance,  which  was  introduced  by  S.  Fan’s  group  at  Stanford  University.  Unlike  usual  guided 
modes,  in  which  light  is  confined  within  the  slab,  guided  resonances  provide  a  unique  way  to 
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channel  light  from  within  the  slab  to  the  external  environment.  Conversely,  light  can  be  coupled 
into  the  device  through  the  PhC  as  well.  Furthermore,  the  sharp  Fano-like  guided  resonance  peaks 
offer  a  simple  yet  efficient  way  to  implement  wavelength  selective  devices,  such  as  sensitive 
filters,  which  are  key  elements  in  satellite,  mobile,  wireless,  microwave,  and  millimeter-wave 
communication  systems.  In  addition,  if  two  PhC  slabs  are  combined,  they  form  a  sensitive 
displacement-tunable  structure  that  can  be  utilized  in  designing  functional  tunable  devices,  such 
as  switches,  filters,  modulators,  and  sensors.  The  key  feature  of  this  tunable  device  is  its 
compaction  and  sensitivity  due  to  the  small  mechanical  motion  required  for  device  operation. 


Frequency  (GHz) 


(a)  (b) 

Fig.  27:  (a)  Schematic  picture  of  a  displacement  sensitive  PhC  filter,  (b)  Experimentally  measured 
transmission  spectra  of  the  GRF;  dashed  and  solid  lines  represent  the  measured  and  3D-FDTD 
calculated  transmission  spectra,  with  the  red  and  the  green  curves  represent  the  gap  as  3mm ,  blue  and 
purple  represent  the  l mm  gap,  respectively. 

However,  most  of  the  studies  on  the  guided  resonance  are  concentrated  on  theoretical 
simulations.  We  demonstrated  a  displacement  tunable  filter,  shown  in  Fig.  27(a).  We  simulated 
and  experimentally  measured  the  transmission  spectra  of  this  structure  using  the  3D-FDTD 
method  and  the  same  measurement  setup  for  measuring  single-slab  based  guided  resonant  filters. 
The  slabs  are  the  same  as  in  the  single-slab  case,  in  which  the  PhC  consists  of  a  square  lattice  of 
air  holes  with  a  -  3.3  mm  and  r  =  1  mm.  The  measurement  is  conducted  after  the  distance 
between  the  slabs  is  changed  from  1  mm  to  3  mm,  which  is  in  the  sub-wavelength  range.  Fig. 
27(b)  shows  the  measured  and  the  3D-FDTD  calculated  transmission  spectra  of  this  tunable  filter; 
the  dashed  lines  represent  the  measured  spectra  and  the  solid  lines  display  the  3D-FDTD 
calculated  results.  The  red  and  green  curves  representing  the  transmission  spectra  with  the 
distance  between  the  two  slabs  equals  to  1  mm,  whereas  the  purple  and  blue  curves  representing 
the  spectra  with  the  distance  between  the  slabs  is  increased  to  3  mm.  From  these  spectra,  one  can 
observe  that  the  wave  totally  transmitted  at  3  mm  and  highly  reflected  at  1  mm.  Furthermore,  the 
30dB  drop  rate  is  very  promising  for  various  applications,  including  sensors. 
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3.0  DISCUSSION  AND  CONCLUSION 


To  summarize,  we  designed,  fabricated  PhCs  for  millimeter  applications  and  experimentally 
demonstrated,  for  the  first  time,  self-coil imation  effect  using  both  amplitude  and  phase,  for  low 
index  contrast  photonic  crystals.  Theoretically,  for  high  index  contrast  photonic  crystals  self- 
collimation  is  easier  to  be  achieved  because  it  is  easier  to  engineer  the  EFCs  to  desired  shape  and 
it  is  easier  to  avoid  band  degeneracy.  Consequently,  we  can  conclude  that  self-collimation  is  a 
very  widely  existing  phenomenon  for  PhCs  and  a  new  guiding  and  routing  mechanism  for  various 
applications.  This  work  also  clarified  two  different  guiding  mechanisms  of  PhCs,  namely 
“bandgap  plus  defects”  and  self-collimation,  because  there  is  no  bandgap  at  all  in  our  PhCs.  Due 
to  fewer  limitations  for  index  contrast  materials,  this  new  guiding  and  routing  mechanism  will 
play  a  more  important  role  in  optical  integration  and  millimeter-wave  communications  as  we 
expected.  In  millimeter-wave  regime,  the  possible  development  of  narrow  beam  maintenance  for 
acceptable  size  would  enable  better  tracking,  detection  and  surveillance  in  modem  radars  and 
antennas. 

In  addition,  we  proposed  the  use  of  the  unique  dispersion  property  of  photonic  crystals, 
namely  self-collimation,  to  supply  a  good  feeding  wave  for  a  dielectric-grating  antenna.  With  the 
integration  of  such  a  photonic  crystal,  we  designed  and  fabricated  a  millimeter-wave  grating 
antenna.  In  the  experiment,  we  characterized  the  device  by  observing  the  near  field  distribution  of 
the  integrated  device.  Experimental  results  showed  that  self-collimation  photonic  crystal  is  a  good 
candidate  for  feeding  the  grating  antenna.  Both  the  photonic  crystal  and  the  grating  antenna  have 
compact  dimensions,  which  make  them  easily  integrated  with  other  planar  devices.  In  many 
applications  the  use  of  navigation  and  communication  waveguides  and  antennas  renders  an 
otherwise  structurally  sound  platform  at  nsk  of  being  compromised.  However,  the  antenna  system 
we  developed  preserves  the  structural  integrity  of  naval  platforms  as  a  well  as  maintain  low 
scattering  cross-sections. 

Furthermore,  a  3D  simple  cubic  PhC  was  fabricated  and  its  dispersion  properties  were 
investigated.  Experiments  demonstrated  3D  self-collimation  in  the  PhC  by  mapping  3D  field 
distribution  inside  the  PhC.  This  is  the  first  time  that  a  “non-diffractive  material”  based  on  3D 
PhC  has  been  fabricated  and  characterized.  Therefore,  these  results  have  both  theoretical  and 
application  significance,  and  constitute  a  breakthrough  toward  the  development  of  novel 
electromagnetic  materials. 

Finally,  we  experimentally  validated  hybrid  lattice  waveguides,  channel  izers,  and 
displacement  tunable  filters.  The  demonstration  of  the  hybrid  lattice  wavelength  selective  devices 
indicated  that  this  hybrid  lattice  structure  is  not  only  a  potential  candidate  for  filter  and 
channelizer  applications,  but  also  opens  a  new  direction  for  designing  devices  by  utilizing  a 
combination  of  different  properties  and  lattices  of  PhCs.  In  addition,  we  presented  displacement 
tunable  filters  for  filter,  sensor,  and  switch  applications.  The  demonstration  of  channelizers  and 
displacement  filters  suggested  that  these  ultra-small  size  devices  have  great  potential  in 
developing  millimeter-wave  signal  processing  systems  with  wide  bandwidth  and  high  processing 
capabilities.  These  results  offer  opportunities  in  designing  sensors  with  the  directly  signal 
processing  abilities.  To  this  end,  we  will  continue  our  work  on  developing  high-speed,  compact 
and  low  cost  millimeter-wave  signal  processing  devices  and  systems. 
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14:  (a)  The  equi-frequency  surfaces  for  frequency  f=  1 2.3GHz.  The  third  and  fourth  bands  are 
degenerate,  (b)  The  equi-frequency  surface  of  the  fourth  band,  (c)  The  equi-frequency  surfaces  shown 
as  a  series  of  contours  for  different  kz  components.  The  black  circle  on  each  diagram  is  the  dispersion 

contour  for  air . 18 

15:  (a)  The  FWHM  beam  profde  in  the  PhC  at  f—  1 2.2GHz.  (b)  The  evolution  of  beam  diameter  in  the 

PhC  and  common  materials . 20 

16:  (a)  The  relationship  between  the  attenuation  and  working  frequency,  (b)  The  attenuation  along  the 

z-ctxis  at  f=  1 2.4 GHz . 21 

1 7:  The  hybrid  lattice  structure  with  a  rectangular  lattice  of  air  holes  in  the  middle  surrounded  by  a 

triangular  lattice . 22 

18  :  Dispersion  diagrams  of  the  triangular  and  rectangular  lattices . 23 

19:  EFCs  of  the  rectangular  lattice:  (a)  first  band  and  (b)  second  band . 23 

20:  The  2D-FDTD  simulation  of  light  propagating  in  the  PhC  with  a  square-like  EFC . 24 

21:  The  2D- FDTD  simulation  of  light  propagating  in  the  PhC  with  a  square-like  EFC . 24 

22:  The  schematic  picture  of  the  designed  hybrid  lattice  WDM,  where  a  rectangular  lattice  serves  as 
the  channel  waveguide  and  triangular  lattices  serve  as  the  boundary.  Point  defects  are  introduced  in 

the  triangular  lattice  region  as  the  cavities . . . 25 

23:  (a)  and  (c),  experimental  measured  surface-scattered  fields  at  12  and  10.8  GHz;  (b)  and  (d),  2D- 
FDTD  simulations  of  steady-state  Hz  fields  at  12  and  10.8  GHz;  (e),  fabricated  HLPWDM;  (f),  2D- 
FDTD  calculated  transmission  spectrum  at  the  cavity . 26 
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Fig.  24:  2D-PWM  calculation  of  the  dispersion  diagram  of  the  hybrid  lattice ,  structure  field  patterns,  and 

phase  at  the  10.8  and  the  12  GHz . 27 

Fig.  25:  Schematic  of  the  J  x  4  WDM  based  on  the  hybrid  lattice . 27 

Fig.  26:  2D-FDTD  calculated  steady  state  Hz  fields  of  1  x  4  hybrid  lattice  PhC  WDM  with  the  10.8 

GHz,  l  LI  GHz,  1 1.4  GHz,  and  12.2  GHz  wavelengths  of  light  dropped  at  different  channels . 28 

Fig.  27:  (a)  Schematic  picture  of  a  displacement  sensitive  PhC  filter,  (b)  Experimentally  measured 
transmission  spectra  of  the  GRF;  dashed  and  solid  lines  represent  the  measured  and  3D-FDTD 
calculated  transmission  spectra,  with  the  red  and  the  green  curves  represent  the  gap  as  3mm ,  blue  and 
purple  represent  the  l  mm  gap,  respectively . 29 
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